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Abstract 

We present a systematical study on the kaon electromagnetic form factors F K ± K o ko(Q 2 ) within 
the kj- factorization formalism, where the transverse momentum effects, the contributions from 
the different helicity components and different twist structures of the kaon light-cone (LC) wave 
function are carefully analyzed for giving a well understanding of the hard contributions at the en- 
ergy region where pQCD is applicable. The right power behavior of the hard contribution from the 
higher helicity components and from the higher twist structures can be obtained by keeping the kx 
dependence in the hard amplitude. Our results show that the kx dependence in LC wave function 
affects the hard and soft contributions substantially and the power-suppressed terms (twist-3 and 
higher helicity components) make an important contribution below Q 2 ~ several GeV 2 although 
they drop fast as Q 2 increasing. The parameters of the proposed model wave function can be 
fixed by the first two moments of its distribution amplitude and other conditions. By varying 
the first two moments af (lGeV) and cJ£(lGeV) with the region of 0.05 ± 0.02 and 0.10 ± 0.05 
respectively, we find that the uncertainty of the kaon electromagnetic form factor is rather small. 
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I. INTRODUCTION 



The electromagnetic form factors provide useful information concerning the internal struc- 
tures of the mesons, and they also provide useful platforms to check the Tightness of the 
perturbative QCD (pQCD) theory. Recently, the electromagnetic form factor F W (Q 2 ) has 
been restudied in Refs.[l|, |2|, yj. It was shown that when all the power suppressed contribu- 
tions, which include higher order in a s , higher helicities and higher twists in the light-cone 
(LC) wave function, and etc., have been taken into account, then the hard contributions 
can fit the present experimental data well at the energy region where pQCD is applicable. 
By comparison the behavior of the kaon electromagnetic form factor is less certain both 
experimentally and theoretically. The kaon is composed by two quarks with different quark 
masses, therefore it becomes a little more complicated to obtain its LC wave functions and 
to compute its electromagnetic form factor. For example, the kaon electromagnetic form 
factor has been studied in Refs.0, 0] in the light-cone quark model only within the soft 
region. Here we shall present a systematical study on the charged/neutral kaon electromag- 
netic form factors in the intermediate and large energy region within the kx factorization 
formalism by properly taking the SUf(3) breaking effects into account. 

The kaon electromagnetic form factor can be obtained through the definition 

(K(p')\J,\K(p)) = (p + p'),F K (Q 2 ), (1) 

where K stands for K ± , K° and K° respectively, the vector current = J2i ^iQil^Qi, with 
the quark flavor i and the relevant electric charge e«. The momentum transfer q 2 = —Q 2 = 
(p — p') 2 is restricted in the space-like region. In the LC quantization and by using the 
Drell-Yan-West (q + = 0) frame |6|], the kaon electromagnetic form factor can generally be 
expressed as 

F K (Q 2 ) = *® £ = ^^ e i [i dx i]i d ^n%(XiM±Ai)^n(Xi,^X + SiqxAi)^ (2) 
n,\i j 

where the summation extends over all quark/gluon Fock states which have a non- vanishing 
overlap with the kaon, Cj is the electric charge of the struck quark, \I/ n are the corresponding 
wave functions which describe both the low and the high momentum partons, [dxj][dkjj_] n 
is the relativistic measure within the n-particle sector and 5, = (1 — Xi) or (— xi) depending 
on whether i refers to the struck quark or a spectator, respectively. 
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Similar to the pionic case, it can be found that the nominal power law contribution to 
F K (Q 2 ) as Q -y oo is F K (Q 2 ) ~ l/(Q 2 ) n - 1 in the light-cone gauge (A+ = 0) 7^ under 
the condition that n quark or gluon constituents are forced to change direction. Thus only 
the qq component of \I/(( 1 ~ a: ) ( 3)(a; ; kj_, A) contributes at the leading 1/Q 2 . For the large Q 2 
region, the hard contribution to the kaon electromagnetic form factor can be written as 

F K {Q 2 ) = 

]T ej J [dx] [dy] [d 2 k ± ] [dH^-^ix, k ±J X)T H (x, y, <u, k ± , 1 ± , A, A')^ 1 "^, U, A') + ■ (3) 

3 

where the ellipses represent the higher Fock states' contributions, [dx] = dx\d2S(l —x\ —x^) 
and [d 2 kj_] = <i 2 kj_/167r 3 . ^J^ 1 ~ x ^(x,h±, X) is the valence Fock-state LC wave function 
with helicity A and with a cut-off on |kj_| that is of order (1 — x)Q. Such a cut-off on |kj_| is 
necessary to ensure that the wave function is only responsible for the lower momentum re- 
gion. And the hard scattering amplitude T# contains all two-particle irreducible amplitudes 
for 7* + qq — > qq. 

The LC wave function provides useful links between the hadronic phenomena in QCD at 
large distance (non-perturbative) and small distance (perturbative). A LC wave function is 

nn 

a localized stationary solution of the LC schrodinger equation idy&(r)) = H lc \^>{t)) [8|, |9|, 
which describes the evolution of a state |\£(t)) on the LC time r = x + = x° + x 3 in the 

caon, its LC wave functions 

3, 



physical LC gauge A + = 0. For the valence quark state of the 
can be defined in terms of the bilocal operator matrix element 



(K(p)\q p (z)q a \0) 



^75^0, kj.) - fiKl5 



%{x, k±) - n>*n u - p" 



6 6 <9k i 



>(4) 

a/3 



where \lk is the phenomenological parameter: fix = Mj^/(m s + m u ) for K and \ik = 
M\ I (m s + rrid) for K° or K° respectively, which is a scale characterized by the chiral 
perturbation theory. ^/ K (x,'k±) is the leading twist (twist-2) wave function, kj_) and 

ty a (x, kj_) are sub-leading twist (twist-3) wave functions that correspond to the pseudo- 
scalar structure and the pseudo-tensor structure respectively. The wave function *$?(x, kj_) 
stands for ^* k, and ty a respectively) satisfies the normalization condition 
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where the decay constant fx = 160MeV The distribution amplitude (DA) (j)(x) and 
the wave function ^(x, kj_) are related by 

, . , 2y/6 r d 2 k± . . 

f K J|k ± |<M/ 167T 3 

Non-leading perturbative contributions to the kaon electromagnetic form factor include 
the higher order in a s , higher helicities and higher twists in the LC wave function, and etc. 
Similar to the pionic case l|, [2] , it is substantial to take kr dependence in the wave function 
into account and to keep the transverse momentum dependence fully in the hard scattering 
amplitude in the kx factorization formalism within the LC framework. In present paper, 
we shall calculate all the helicity components' contributions to the kaon electromagnetic 
form factor within the LC pQCD framework, which is consistent with the using of LC wave 
function. Another important power correction is from the higher twist structures in the kaon 
DA. The end-point singularity becomes more serious for the higher twist structures, because 
the asymptotic behavior of the twist-3 DAs, especially <f)p S {x) = 1, so the calculations for 
these higher twist contributions have more uncertainty than that for the leading twist. It 
means that one should use the twist-3 wave function with a better behavior in the end-point 
region than that of the asymptotic one so as to give a more reliable estimation of the higher 
twist structures' contribution. The Brodsky-Huang-Lepage (BHL) prescription provides a 
useful way to construct a wave function with better end-point behavior |12j , we shall adopt 
it to construct the kaon LC wave functions for the present purpose, and then we discuss its 
uncertainty for the kaon electromagnetic form factor. The 5'?7/(3)-breaking effects shall also 
be included for constructing the kaon LC wave function. 

The reminder of the paper is organized as follows. Sec. II is devoted to present the main 
properties of the kaon electromagnetic form factor and the formulae for the twist-2 and twist- 
3 contributions to the kaon electromagnetic form factor within the kx factorization approach. 
Numerical results for the kaon electromagnetic form factor are presented in Sec. III. The last 
section is reserved for conclusion and summary. 
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FIG. 1: Six leading order time-ordered Feynman diagrams for the hard scattering amplitude Th, 
where p 1 = (a;i,kj_), p 2 = {x 2 ,-k ± ), p[ = (yi,yiq± + p' 2 = {y2,V2<i_L - U)- 

II. KAON ELECTROMAGNETIC FORM FACTOR WITHIN THE k T FACTOR- 
IZATION FORMALISM 

Because, K + = us and K~ = su, one can find that F K + = —F K - according to Eq.(j2]). 
Similarly, since K° = ds and K° = sd, it can be found that F K o = —F^o. So, we only 
need to calculate the K + and K° form factors, where e\ = 2/3, e 2 = 1/3, mi = m u and 
m 2 = t^s for F K + and t\ = —1/3, e 2 = 1/3, mi = and m 2 = m s for F K o respectively. 
Here m u = ^ m s stand for the light constitute quark masses. Further more, in doing the 
calculation of the hard scattering amplitude with the fcy factorization formulism, we shall 
treat the current quark masses of u, d and s to be zero due to their smallness in comparison 
to the involved hard scale. Then the calculation procedure for the hard scattering amplitude 
is the same as that of the pionic case 



A. Formulae for the twist-2 contribution to the kaon electromagnetic form factor 

In the intermediate and large energy region, one can apply the pQCD approach and use 
the valence Fock state to estimate the kaon electromagnetic form factor since the applica- 
bility of pQCD in the intermediate and large energy region has been proved by involving 



the kx dependence [13| . The lowest-order contribution for the hard scattering amplitude T# 
comes from the one-gluon exchange Feynman diagrams as shown in Fig.([T]). To simplicity 
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TABLE I: Full form of the LC wave function ^(x,h±, X) = (p[x,\c±)x- ^f(x, kj_, A) stands for 
VI/rt (x,kj_, A), ty p (x, kj_, A) and ^ a {x, kj_, A), respectively. 



A1A2 


TT 


TI 


^AiA 2 (^,k_L 


,A) 


(a 1 +a 2 )(k x -iky) , 


k±) 


(aia 2 -k^) / , x 


[2(a?+k* )(^+k^ (IV^^l- 1 


[2(^+k^)(a|+k^)]l/2^^' k ^^ 


A1A2 


IT 


II 


^ Al A 2 (a;,k_L 


,A) 


(aia2-k^) ^ 




(ai+a 2 )(fca : +iA:a) / , n 


[2(a^+ki )(a^+ki )]Va^^ 


[2(aJ+ki)(a|+ki)]V2V^' k± '' 



our notations, we separate the spin-space wave function x \ x -> k±, A) out from the whole LC 
wave function, i.e., ^^~ x ^ Q (x, kj_, A) — > kj_, A)<£>( 1_a: ^(x, kj_), where vE^ 1-2 ^^, kj_) 

is the light-cone wave function of the valence Fock state with a cut-off |kj_| of order (1 — x)Q 
and the spin space wave function x K ( x > kj_, A) that comes from the spin space Wigner ro- 



tation can be found in Ref. 14| 1 , which is given in TAB]!] One can combine the spin-space 



wave function x \ x i kj_, A) into the original Tjj to form a new one, i.e., 

T H = (er6 + e 2 6)7f 1+A2=0) (TI-TI) + (e x 6 + e 2 6)Ti, Al+A2=0) (IT-IT) + 

( ei £ + e 2 e 2 )4 Al+A2=1) (TT-TT) + (eitf + e2^)4 Al+A2=_1) (II-II) , (7) 

where Ai )2 are the helicities for the (initial or final) kaon's two constitute quarks respectively. 
It is found that there is no hard scattering amplitude with quark and antiquark helicities 
being changed due to the fact that the quark helicity is conserved at each quark-gluon 
(photon)-quark vertex in the limit of vanishing quark mass, is the electric charge of the 
struck quark, and £• are coefficients derived from x K { x i kj_, A), 

(aia 2 - k5_)(aia 2 - 



6 2[(a? + ki)(a| + ki)(af + li)(a 2 2 + li)]V»' 

(ai + a 2 )(a'i + Q 2 )(k_L ■ 1± + ikj_ x 1 ± ) 
^ := 2[(a? + ki)(a| + ki)(af + li)(a 2 2 + li)]V2' 

P (^-kD^-ij) 

v ' " 2[(6? + ki)(6| + ki)(6f + li)(6' 2 2 + li)] 1 / 2 



and 



(&i + 62)(6 / i + & , 2 )(k ± -U + »k ± xU) 
^ 2 2[(fc 2 + ki)(6 2 + ki)(6f + li)(6 2 2 + li)]V2 ' ^ 



1 Setting mi = to 2 = m q , we return to the results of the case of pion. 
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where 



ai = xM a + mi, a 2 = (1 — x)M a + m 2 
a[ = yM' a + mi, a' 2 = (1 - y)M^ + m 2 
61 = xM 6 + m 2 , b 2 — (1 - x)M b + m x 
6; = yM^ + m 2 , 6' 2 = (1 - y)M' b + m x 



M 2 = ^ 



/// 



x 



I + k l 
1 -x 



M' 2 = 



m; 



y 



+ 



m 



M 6 2 = — H } -, M b = — + 1 



x 1 — x y 1 

Consequently, the above coefficients can be further simplified as 

6 



and 



£2 



[mi(l 


— x) + m 2 x] [mi ( 1 — y) + m 2 y] 


2 v /A; 2 +[m 1 (l 


- x) + rn 2 x\ 2 \Jt]_ + [mi(l - y) + m 2 y] 2 




(k ± • 1 ± + ik ± x 1 ± ) 


2^1 + ^(1 


- x) + m 2 x} 2 \Jl 2 ± + [mi(l - y) + m 2 y] 2 


[m 2 (l 


- x) + mix] [m 2 ( 1 - y ) + m x y] 


2^*1 + Mi 


- x) + m x xf\jl\ + [m 2 (l - y) + miy] 2 




(k ± -l ± + ik ± x 1 ± ) 



2 \Ao. + [^2(1 - x) + m x xf\jl\ + [m 2 (l - y) + m x y] 2 
Schematically, the total hard scattering amplitude can be written as 

T H = (ei6 + e 2 6)4 Al+A2=0) + [ei(6 + &) + e 2 (£ 2 + £ 2 *)]4 Al+A2=±1) 



with 



,(Ai+a 2 =o) . lQ7rC F a s (fi 2 j 



(10) 



(11) 



H 



x (((x - l)q ± 2 - 2k ± • qj(21 ± • q ± + (y - l^ 2 ))" 1 



(1 - x)(l - y)xy 
((x - 1)(21_l • q ± + (y - l)q ± 2 ) - 2(y - l)k ± • q^ 1 x 

2(y - 1)3/(1 - y + x(2y - l))(k ± • q ± ) 2 + (x - l)x(21 ± • q ± + (y - l)q ± 2 ) 
((1 - y + x(2y - 1))(1 ± • qj + 2(x - l)(y - l)yq ± 2 ) - 



(x - l)(y - l)y(k ± • q ± ) • (8x(l ± • qj + (1 - y + x(6y - 5))q ± 2 ) 



(12) 



and 



n(Ai+A 2 =±l) 
if 



87rCVa s (/i 2 ) 
(1 - x)(l - y)xy 



x (((x - l)q ± 2 - 2k ± • q ± )(21 ± • q ± + (y - l^ 2 ))" 1 
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{(x - 1)(21jl -<u + (y- l)qx 2 ) - 2(y - l)k ± • qx)" 1 



2(x - l)a;(l ± • q ± ) 



(y - l)(2y(k ± • q ± ) 2 + (x - l)(x(l ± • qj - y(k ± • q ± ))q ± 5 



(13) 



where the scale /i? = Q 2 . It can be found that the leading contribution from the higher 
helicity components is of order l/Q 4 , which is next-to-leading contribution compared to that 
of the ordinary helicity components. 

With the help of Eq.Q, we can obtain the leading-twist hard part contribution to the 
kaon form factor. And after integrating over the azimuth angles for k^ and we obtain 
the contribution from the usual helicity components (X\ + A 2 = 0), 



rfu>ist2,(Ai+A 2 =0) ^Q2\ 



(ei^i)C F a s (//J)|k ± ||l ± | 
dxdyd^d^ oo — J ip(x, k ± )ip [y,l±) x 



32-K s xy 



x(x + y — 1 — 2xy) y(x + y — 1 — 2xy) 



x + y — x 2 — y 2 



+ 1 ei <-> e 2 , m 1 <-> m 2 



(1 - y)y/l- V % (1 ~x)(l- y^l-rftjl-rft 



(14) 



and the contribution from the higher helicity components (Aj + A 2 = ±1) 



itwist2,(\ 1 +\ 2 =±l) (r .2 



K 



dxdydr\\dT]\ 



(e 1 ^ 3 )C , Fa a (/i?)|kj.||lj 



QAir 3 xy 



-<p(x,k ± )<p*(y,l ± ) 



x 



(x + y — 2xy)- 



l-Jl^vl 1- yfl^vl 



1 - x)(l - y)r)ir} 2 Jl - r/f a/1 - iff 



+ < ei <-> e 2 , mi « m 2 , 



(15) 



where without loss of generality, we have implicitly assumed that the radial kaon wave func- 
tion <p(x, kj_) depending on kj_ through k\ only, i.e. <p(x, kj_) = <p(x, k\). The terms in the 
big brace are obtained by transforming the terms out of the brace through the transforma- 



IjLkHkl 



f , |kj_| = Q(l-a0t7i/2 



tion ei <-> e 2 and mi <-> m 2 . £3 = , = — ,— - - -- ■ = 

S y/kl + [rn 1 (l-x)+rn 2 x} 2 y/ll + [m 1 (l-y)+m 2 y} 2 

and = Q(l — y)7] 2 /2, with r/ 12 in the range of (0, 1). An overall minus sign in Eq. (|15p 
implies that the higher helicity components shall always suppress the contribution from the 
usual helicity components. 
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B. Formulae for the twist-3 contributions to the kaon electromagnetic form factor 



The end-point singularity becomes more serious for the higher twist structures, because 
the asymptotic behavior of the twist-3 DAs, especially <fip S (x) = 1, so the calculations for 
these higher twist contributions have more uncertainty than that for the leading twist. As 
has been pointed out in Ref . [ijj] , after including the parton transverse momenta, large dou- 
ble logarithmic corrections a s In 2 k± and a s In 2 x appear in higher order radiative corrections 
and can be summed up to all orders. The relevant Sudakov form factors from both k± and 
the threshold resummation can cure the endpoint singularity and then the main contri- 
bution comes from the perturbative regions. For the present purpose, it is convenient to 
transform the kaon form factor into the compact parameter b space. In the large Q 2 region, 
by considering only the lowest valence quark state of the kaon and by doing the Fourier 
transformation of the wave function with the formula, 

r d 2 b 

tfCa^kj.;^) = J ^e- ,b - kl $(ij,b;/i / ) 1 
we can transform the kaon electromagnetic form factor into the compact parameter b space, 



F K (Q 2 ) = I [rfx i (ib][ci|/ J dh]^(x i ,b;/i / )T(a; i ,b;y i ,h;/i / )^(?/ j ,h; y u / ) x S t (xi)S t (yj) x 

exp(-S(xi,yj,Q,b, h; ///)), (16) 



where [if = ln(fj,f/ Aqcd), [dxidb] = dxidx2d 2 b5(l — x\ — X2)/(167r 3 ) and the hard kernel 

f(x i ,b; % -,h;/x / ) = J |^ ^ 2 e- ih ^-^T(x u k ±i ; Vj , l ±j ; [i f ). 

The factor exp(— S(xi, jjj, Q, b, h; [if)) contains the Sudakov logarithmic corrections and the 
renormalization group evolution effects of both the wave functions and the hard scattering 
amplitude, 



S(x 1 ,y 1 ,Q,h,h;[i f ) 



(E s ^ b,Q)+£, h, Q) ) - i In U - 1 In h 

\i=i j=i J Pi b Pi h 



(17) 



where b = ln(l/6AQc , r>), h = \n(l/hAQCD) an d s(x,b,Q) is the Sudakov exponent factor, 
whose explicit form up to next-to-leading log approximation can be found in Ref. jig . In 
Eq. fflBT) . S t (xi) and S t (yi) come from the threshold resummation effects and the exact form of 



each involves one parameter integration 17j . In order to simplify the numerical calculations, 
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we take a simple parametrization proposed in Ref.[l 

2 1+2 T(3/2 + c 



S t (x) 



[x(l-x)Y 



(18) 



v^r(i + c) 

where the parameter c is determined around 0.3. 

With the help of the above equations, we obtain the formula for the twist-3 contributions 
to the kaon electromagnetic form factor, 

1287T/4 n 



rptwisfi 



(Q 



dxdy / bdbhdha s (jjf)T(x } b;y,h.;/j,f)St(xi)St(yj 



x 



Vp(x,b;n f )- 



exp [S(x i ,y j ,Q,b, h; fi f )\ 



(19) 



where x = (1 — x), y = (1 — y) and i>'*(y, h; [if) = d^*(y, h; fif)/dy. The hard scattering 
amplitude T(x, b; y, h; fif) is given by 

f(x, b; y, h; N ) = ^ (v^Q&) (*(& " W A) (vW) + 

- b)K (VVQh) h (VVQb) ), (20) 

where the higher power suppressed terms such as (k± 2 /Q 2 ) has been neglected in the nu- 
merator, 1$ and Kq are the modified Bessel functions of the first kind and the second kind 
respectively. To ensure that the pQCD approach is really applicable, one has to specify 
carefully the renormalization scale iif in the strong coupling constant. Here we take the 
scheme that is proposed in Refs. 18j, i.e. its value is taken as the largest renormalization 
scale associated with the exchanged virtual gluon in the longitudinal and transverse degrees, 

/if = max(y/x r §Q, 1/6, 1/h). (21) 

The full form of the kaon LC wave function have four helicity components (Table. [T]): 
namely, 

* = (* n ,* n ,* it ,* u ), = * ff ) (22) 

By including the higher helicity components, Eq. (1191) can be improved as 
128vr/4 n 



^twist3 
K 



(Q 



dxdy / bdbhdha s (/jf) x T(x } b;y,h; /j,f)S t (xi)St(yj 



x 



y V(%, Ai, A 2 ) + E V(*' a , Ai, A 2 ) + i E A a ) 

A1A2 A1A2 A1A2 

x exp [-S^a*, y,-, Q, b, h; ///)] , 



(23) 
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where ^f' a — d^f a /dx and 

E V(%, A 1; A 2 ) = (*; n * pn + *;; T *pir) - (*;n*PTT + Ku^pu)' 
A1A2 

E ^(^, Ai, a 2 ) = (#; n *; n + *; iT *; iT ) - (*; n *; n + *; u *' ctU ), 
A1A2 

A1A2 

For the hard scattering amplitude T(x, b; y, h; ///), we have implicitly adopted the approx- 
imate relation, i.e. T(x, b; y, h; ~ — T(x, b; ?/, h; yU/)^ +J - T , since it can be found 
that 

f(x, b; y, h; ^)^+U = _f (a;, b; h; ^)^t + 0(1/Q 2 ). (24) 

III. NUMERICAL RESULTS FOR THE KAON ELECTROMAGNETIC FORM 
FACTOR 

Based on the formulae presented in the last section, we discuss sequentially the leading 
and the power suppressed contributions to the kaon electromagnetic form factor within the 
space-like region in the following. The differences between the pion and kaon electromagnetic 
form factors shall also be discussed in the due places. In the numerical calculations, we use 
A^ -4 '' = 250MeV. As for the phenomenological parameter /j,k, which is a scale character- 
ized by the chiral perturbation theory, we take its value to be \xk — 1.70 GeV. And for def- 
initeness, we take the conventional values for the constitute quark masses: m u ^ = 0.30GeV 
and m s = 0.45GeV. 

A. LC wave function of the kaon 

In order to obtain the numerical results for the kaon electromagnetic form factor, we 
need to know its LC wave functions. One useful way is to use the approximate bound state 
solution of a hadron in terms of the quark model as the starting point for modeling the 
hadronic valence wave function. In combination of the spin-space wave function y that 
comes from the Wigner rotation pjj, the full form of the kaon LC wave function can be 
written as, ty(x,k±, X) = (p(x,~k±)x- The explicit form of the spin-space wave function \ 
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can be found in TABJU As for the radial part of the wave function, we adopt the model 



constructed in Refs.[2( 



21(, which is based on the BHL-prescription 
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1 + B K C 3 1 /2 (2x - 1) + C K Ct'\2x - 1)] 



3/2, 



A 



K 



x(l — x) 



<p a (x,k±) 
(p p (x, ki_) 



x exp 
A a exp 



-P 



K 



-Pi 



[1 + B p C\' 2 {2x - 1) + C p C 2 '\2x - 1)] 



k\ + m\ k\ + m 2 

x 1 — x 

k\ + m\ kj_ + m{ 



+ 



x 



1 

1/2, 



X 



(25) 
(26) 



A, 



x(l — x) 



x exp 



-P. 



k 2 , + m? /c 2 + m 2 , 



A" 



+ 



(27) 



a; 1 — x 

m i — m d and m 2 = m s for F K o(Q 2 ). C^ 2 (2x— 1) 



where mi = m„ and m 2 = m s for Fk+{Q 2 ), 
and C*y 2 (2x— 1) are Gegenbauer polynomials. A more complicated model that is also based 
on BHL-prescription is suggested in Ref.j^]. Numerically it can be found that the two model 
wave functions behave very likely under the same constraints. Additionally, it has argued 
in Ref. 



22j that an extra factor Jdk z /dx should be included into the LC wave function, 



otherwise, one can not obtain the right asymptotic behavior for the pion electromagnetic 
form factor at Q 2 — > oo. However, we have checked that without this factor, one can still 



obtain the right power behavior for the pion electromagnetic form factor as shown in Ref. 



a, 



and such a factor will not bring any new features for the LC wave function if its parameters 

20, Q to 



are determined properly. So we shall adopt the simpler form as suggested in Refs. 
do our calculation. 

The four parameters Ak, Bk, Ck and /3k of <Pk(x, kj_) can be determined by the first 
two Gegenbauer moments af and af- of 4>k{x), the constraint (k^_)]^ 2 w 0.350GeV [23] and 
the normalization condition 

OiCt2 



1 , f d 2 k ± a x a 2 - k 2 ± 

o JkK,l 16vr3 [(a 2 + ki)(ai + ki)]V2^ x ' k ^ 

1 , f d 2 k ± m 1 (l-x)+m 2 x . , . 

dx 2 2 1^3 r— = </W, k ±) 

Jk\<iil 167T' 3 



(28) 



1 k^ + [mi(l ~ x) + m 2 xY 2v6 
where /io stands for some hadronic scale that is of order 0(lGeV). Here, the wave function 
is normalized to /ft-/2v / 6 only for convenience, which is different from that of Refs. 2(J, [21] 
that is normalized to one 2 , where the factor f K /2\/6 has been absorbed into the hard part 



It is noted that the unit of Ak in Ref.[20( should be corrected from GeV 1 to GeV 
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of the B — > K transition form factor. And the average value of the transverse momentum 
square of kaon is defined as 



/ dxd 2 h± 


kil 


\^k(x, k ± )| 2 / dxd 2 ~k± 


kil 


(/?jc(x,k_L)| 2 


J dxd 2 \t± 


^(X,k ± )|2 fdxd*k ± \<p K (x,k ± )\* 



The first Gegenbauer moment has been studied by the light-front quark model 24j, the 



LCSR approach 



25 



26 



27 



28 



29] and the lattice calculation 30j, |3l|, ]32) and etc. 



higher Ge gen bauer moments, such still determined with large uncertainty 



27, 



33 



34 



2a 



he 



351 ] . In the following calculation, if not specially stated, we take af (lGeV) = 0.05 



n 



26j and af- (lGreV) = 0.115 |35j to be their default values, and shall discuss the uncertainties 
caused by these two factors in due places. By taking these default values, we obtain 

A K = 12.55GW 1 , B K = 0.0605, C K = 0.0348, (3 K = 0.8706GW" 1 . (29) 

The parameter A a of tp a {x, kj_) can be determined by its normalization condition similar 
to Eq. (j28l) . i.e. A a = 65.04Gel / ~ 1 . And the coefficients A p , B p and C p of (p p (x,\t±) can be 
determined by the DA moments of <p p (x, kj_). To discuss the uncertainty caused by ty p , we 



take two groups of DA moments that have been obtained in Refs 



25 



coefficients A p , B p and C p , where the moments in Ref.[25| are derived by using the QCD 



light-cone sum rules and the moments in Ref. 
action from the instanton: 



36j to determine the 



36j are derived based on the effective chiral 



Group 1 [25] : (x )* = 1, (a; 1 )^ = 0.06124, (x 2 ) 



h 



Group 2 [36] : (x )* = 1, (x 1 )* = 0.00678, (x 2 ) 



1\K 



„2\K 
fp 



0.36757, 
0.35162. 



(30) 
(31) 



Here the moments are defined as (x l ) p = f dx{2x_ — 



25 



)V P (1 - x,fi ) with i = (0,1,2). It 



x, (jl ) other than 



361 ] are for 



should be noted that the moments defined in Ref. 
(t> p (x,Ho), since in these references x stands for the momentum fraction of s-quark in the 
kaon (K), while in the present paper x stands for the momentum fraction of the light quark 
q in the kaon (K). Taking the above two groups of DA moments, the parameters of <p p (x, k±) 
can be determined as, 



Group 1 : 
Group 2 : 



A 1 = 12.12GeV-\ B\ = 0.3062, Cl 



A 2 = 12.04GeV 



-i 



B 2 = 0.4711, Cl 



1.604, 
1.506. 



(32) 
(33) 



It is found that both distribution amplitudes are double humped curves and are highly sup- 
pressed in the endpoint region. Such feature is necessary to suppress the endpoint singularity 
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coming from the hard-scattering kernel and then to derive a more reasonable results for the 
twist-3 contributions to the kaon electromagnetic form factor. 



B. Valence Fock state contribution in the low energy region 



At the present, the experimental data on the kaon electromagnetic form factor are con- 
centrated in the low energy region Q 2 < lGeV 2 , c.f. Ref. The soft part contribution 
can be written as 



Fk(Q 2 ) = £5> £ dx J ^v^( x , k± ,A)vM*,k' ± ,A') 



(34) 



where A, A' are the helicities of the wave function respectively, and the first term is the 
lowest order contribution from the minimal Fock state (valence Fock state) and the ellipses 
represent those from higher Fock states, which are down by powers of 1/Q 2 and by powers of 
a s in the large Q 2 region. In general, the kaon electromagnetic form factor should sum over 
all of higher Fock state contributions in the low energy region. If only taking the leading- 
twist LC wave function of the valence Fock state, we can examine the contribution from the 
valence Fock state in the low energy region, i.e. 



Ff\Q 2 )=e 1 



dx 



d 2 k. 



■K(p K (x, k ± )(p* K (x,k'_ 



+ e 2 



(35) 



167T 3 ' 

where e\ = 2/3, e 2 = 1/3, m 1 = m u and m 2 = m s for F%1' and e\ = —1/3, e 2 = 1/3, 
mi = rrid and m 2 = m s for FgjP respectively. The terms in the second bracket is obtained 
by transforming the terms in the first bracket through mi <-> m.2, the coefficient k that is 
from the spin-space Winger rotation can be written as 

(aid 2 - k^_)(a / 1 a 2 - k'^_) + (oi + a 2 )(a' 1 + a 2 )k_L • k'j_ 
K ~ [(a? + ki)(al + ki)(a? + k'i)(a' 2 2 + k'i)]V* 

_ [mi(l - x) + m 2 x} 2 + k ± ■ k^ 

\Jk 2 _ + [mi(l — x) + m 2 x] 2 y / k 12 + [mi(l — x) + m 2 x] 2 
where k'j_ = k^ + (1 — a;)qj_ for the final state LC wave function when taking the Drell- 
Yan-West assignment 6]. Since mi ^ m 2 , we have F S ^\Q 2 ) ^ 0, which is different from 
the pionic case, i.e. F^ V \Q 2 ) = because 7r° has equal quark masses. 

We proceed to integrate the transverse momentum k^ in Eq. (1351) with the help of the 

Schwinger a— representation method, 

1 1 



a k r(/«) Jo 



a K - l e~ aA da 



(37) 
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Doing the integration over k^, we obtain 

"1 rd 2 k. 



dx 



o 



16/T 3 



K(p K (x, k ± )<p* K (x,k'j 



A\{x) 



dx I dX- 
o Jo 128tt 2 (1 + A) 3 



Q 2 (x - l)f3 2 K X 2 \ /4(1 - x)x(l + X) 



-Q 2 (l - x)\2 + A(4 + A)) + 8m 2 (l + A) 



4x(l + A) J \ P 2 K 

Q 2 {x-l)(5 2 K X 2 ^ 



h 



Q\l-x) 2 X 2 



Ax{\ + A) 



x exp 



(3 2 K [8m 2 a (l + A) + 8m 2 A(A + 1) + Q 2 (l - x) 2 {2 + A(4 + A))] 



4(1 - x)x{l + A) 



(38) 



where the short notations Ak(x) 



[1 + B K Cl ,2 {2x - 1) + C K Cl'\2x - 1)] 



3/2, 



x(l—x) 



m„ 



m\(l — x) +m\x, m 2 = [mi(l — x) + rri2x} 2 , and J n (n = 0, 1) is the modified Bessel function 
of the first kind. Substituting Eq. fl38|) into Eq. fl35l) and doing the expansion in the small Q 2 
limit, we obtain 



ei 



A 2 K (x) 



dx / dX- 
o Jo 16vr 2 (l + A) 2 



cxp 



2(3 2 K (m 2 a + m 2 X) 
(1 — x)x 



m 2 b (l + X) 



x(l — x) 



2P 



+e 2 < mi 



m 2 



(39) 



where the term [m 2 (l + A)] in the second square bracket comes from the ordinal helicity 
components, while the remaining terms are from the higher helicity components. As for the 
mean square radius {r\)v, we obtain 



(r 2 K ) 



dFf\Q 



v 



-6 



ei 



x 



dx 
o .'o 

'1 — x)x 



Q 2 =0 

dX 



3A 2 K (x)(3 2 K 



32tt 2 x(1 + A) ; 



exp 



Am 2 ) 



P 2 



+ m 2 b (l + X) 



A 



+ e 2 < m 1 



(1 — x)x 
m 2 >■ 



(l-a;)(2 + 4A + A 2 ) 



(40) 



The result for the soft contribution to the kaon electromagnetic form factor is shown 
in Fig.©, where the solid line is for of (lGeV) = 0.05 and af(lGeV) = 0.115 and the 
shaded band is obtained with of (lGeV) G [0.03,0.07] and a^lGeV) G [0.05,0.10]. The 
valence quark contribution is slightly below the experimental data, which means that there 
are still some space for the higher Fock state contributions. With the help of Eq. (138p . we 
can estimate the probability of finding the valance states in the charged/neutral kaon, e.g. 
we obtain (P uj = 0.901 < 1.0), which shows that higher Fock states and higher twist terms 
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0.02 0.04 0.06 0.08 0.1 0.12 0.14 

Q 2 (GeV 2 ) 



FIG. 2: The valence Fock state contribution to the kaon electromagnetic form factor \F K + (Q 2 )\ 2 
in the low energy region, where the experimental data is taken from Ref.[37|]. The shaded band 
is obtained with of (lGeV) £ [0.03,0.07] and a%(lGeV) £ [0.05,0.10], and the solid line is for 
af (lGeV) = 0.05 and af(lGeV) = 0.115. 

should also be considered to give a full understanding of the form factor at the energy region 
Q 2 — > 0. Such probability can be further divided into two parts: (-Pig 1+A2_ °' ) = 0.562) for the 
usual helicity components and (Pi^ 1+A2_±1 ' ) = 0.339) for the higher helicity states. It shows 
that the higher helicity components have the same importance as that of the usual helicity 
components in the soft energy region. It is noted that the higher helicity components' 
contribution to the kaon electromagnetic form factor has also been studied with the LC 
framework in Ref. 5[], where the probability of the leading Fock state is just normalized to 
one and the experimental data on the mean-square radius of charged/neutral kaon are used 
to determine the wave function parameters. As argued above, this simple treatment maybe 
not right, since then the contribution from the valence state can be enhanced and become 
important inadequately 3 . 

As for the charged and neutral mean square radii (r 2 K+ ) v and (r 2 K0 ) V) we obtain 
{r 2 K ±)y = 0.570 fm and (r 2 K0 )v = —0.0736 fm 2 , which is consistent with the Ref. 38], 
while experimentally we have {r 2 K± ) 1/2 = 0.560 ± 0.031 fm and (r 2 K0 ) = -0.077 ± 0.010 fm 2 
ll| . Further more, we give a simple estimation of the uncertainties caused by the two Gegen- 



As has been pointed out in Ref. [2|, the condition for the pionic case is more serious, where the probabilities 
for the value quark state is only about 74%. 
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FIG. 3: Leading twist contribution to the kaon electromagnetic form factor in the intermediate and 
large energy region, where the left is for Q 2 F K +{Q 2 ) and the right is for Q 2 F K o(Q 2 ). The dotted 
line stands for the contribution from the usual helicity (Ai + A2 = 0) components, the dashed line 
stands for the contribution from the higher helicity (Ai + A2 = ±1) components and the solid line 
is the total hard contribution, which is the combined result for all the helicity components. The 
dash-dot line stands for the usual asymptotic result of charged kaon. 

bauer moments of the kaon twist- 2 wave function, e.g. by taking af(lGeV) = 0.05 ± 0.02 
and of (lGeV) = 0.10 ± 0.05, we obtain the probability of finding the valence Fock state 
in the kaon P uS = 0.9011o;qio, and the uncertainties of radii (r 2 K± ) l y 2 = 0.570to;o28 ^ m an d 
( r Ko)v — ~ 0-0736lo;oi4 fni 2 - ft should be noted that by taking different values for af and 
af , all the undermined parameters of the wave function should be varied accordingly. 

C. leading twist contribution to the kaon electromagnetic form factor 

With the help of the LC wave function Eq.( l25|) . we show the leading twist contribution 
to the kaon electromagnetic form factor in the intermediate and large energy region in 
Fig-©, where the contribution from the usual helicity component or from the higher helicity 
components are considered. It is shown that the higher helicity components always suppress 
the usual helicity components' contributions to the kaon electromagnetic form factor. The 
usual asymptotic result of charged kaon, i.e. Q 2 F K +{Q 2 )\ asy = 87r/| : a s (Q 2 ), is also presented 
in Fig. ([3]) for reference. It can be found that the leading contribution of the hard-scattering 
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FIG. 4: Uncertainties of the leading twist contribution to Q 2 F K +(Q 2 ) caused by af and af , where 
the left diagram is for fixed af(lGeV) = 0.05 with af (IGeV) = 0.03, 0.05 and 0.07, and the right 
diagram is for fixed af (IGeV) = 0.05 with af(lGeV) = 0.05, 0.10 and 0.15 respectively. 

amplitude from the higher helicity components is of order l/Q 4 , which is next-to-leading 
contribution compared to the contribution coming from the ordinary helicity component, but 
it shall give sizable contributions to the kaon electromagnetic form factor in the intermediate 
energy region. The net contribution shows the right power behavior Q 2 Fk+(Q 2 )\q2_ >00 — > 
const. In the present work, we have considered the kj- dependence both in the wave function 
and in the hard scattering amplitude consistently within the LC pQCD approach, then our 
results present a right power behavior for the higher helicity components' contributions. 
Secondly, in contrary to the pionic case, we obtain F K o(Q 2 ) = —F^o(Q 2 ) ^ at Q 2 ^ 0, 
which are rightly caused by the S*[//(3)-breaking effect and are strongly dependent on the 
constitute quark masses. 

We take the charged kaon electromagnetic form factor as a concrete example to show 
the uncertain caused by af(lGeV) and af (IGeV), which are varied within the region of 
[0.03,0.07] and [0.05,0.15] respectively. We draw the charged kaon electromagnetic form 
factor in Fig. (131) ; where the left diagram is for fixed af(lGeV) = 0.05 with af(lGeV) = 0.03, 
0.05 and 0.07, and the right diagram is for fixed af(lGeV) = 0.05 with af (IGeV) = 0.05, 
0.10 and 0.15 respectively. Q 2 Fk+(Q 2 ) decreases with the increment of af. From Fig(jlj), 
it can be found that the uncertainty of the form factor caused by af (IGeV) = 0.05 ± 0.02 
is small, e.g. it is about ±5% for q 2 e [2, 20]GeV 2 . And the uncertainty of the form factor 



18 



0.3-| 1 0.05 



0.04 




0.0- 0.00 

1 1 1 1 1 1 1 Li 1 1 1 1 1 1 

5 10 15 20 5 10 15 20 

Q 2 (GeV 2 ) Q 2 (GeV 2 ) 

FIG. 5: Twist-3 contribution to the kaon electromagnetic form factor Q 2 F K +(Q 2 ) and Q 2 F K o(Q 2 ). 
The dotted line stands for the contribution from the usual helicity (Ai + A2 = 0) components, the 
dashed line stands for the contribution from the higher helicity (Ai + A2 = ±1) components and the 
solid line is the total hard contribution, which is the combined result for all the helicity components. 
As a comparison, the twist- 2 contribution is shown in dash-dot line. 

caused by a 2 Y (lGeV) varying within a bigger region [0.05,0.15] is also small, i.e. which is 
about 4% — 9% for q 2 G [2, 20]GeV 2 . Q 2 F K +(Q 2 ) decreases with the increment of in the 
lower energy region q 2 < 6GeV 2 and increases with the increment of aff in the higher energy 
region q 2 > QGeV 2 . 



D. twist-3 contribution to the kaon electromagnetic form factor 

We show the twist-3 contribution to kaon electromagnetic form factors Q 2 Fx+(Q 2 ) and 
Q 2 Fko(Q 2 ) in Fig. ([5]), which are obtained with the full form of the LC wave functions 
*&p(x, kj_) and *ff£(x, kj_) and with the Group 1 parameters for ^ p (x, kj_). It is found that at 
the twist-3 level, the higher helicity components' contributions to the form factor are negative 
and small in comparison to that of the usual helicity components. The twist-2 contribution 
is also presented in Fig. ([5]) for comparison. At the twist-3 level, both the charged and the 
neutral kaon electromagnetic form factors decrease with the increment Q 2 , and the charged 
form factor becomes smaller than the twist-2 contribution at around Q 2 = 7GeV 2 , which is 
changed to be Q 2 ~ YlGeV 2 for the neutral case. This implies that the twist-3 contributions 
are sizable in the intermediate energy region and are rightly power suppressed to the twist-2 

19 




0.04 




Q 2 (GeV 2 ) 



Q 2 (GeV 2 ) 



FIG. 6: Uncertainty caused by two types of twist-3 wave function to the electromagnetic form 
factors Q 2 F K +{Q 2 ) (Left diagram) and Q 2 F K o(Q 2 ) (Right diagram). The solid line and the dotted 
line are for Group 1 and Group 2 parameters respectively. 





Q 2 (GeV 2 ) 



Q 2 (GeV 2 ) 



FIG. 7: Uncertainty caused by the parameter (3k of ^ p with Group 1 parameters to the elec- 
tromagnetic form factors Q 2 Fx+{Q 2 ) (Left diagram) and Q 2 F K o(Q 2 ) (Right dia gram). The solid 
line, the dashed line and the dotted line are for (3 K = 0.85GeU _1 , 0.87GeU _1 and 0.89GeU _1 
respectively. 

contributions in the large energy region, which is similar to the pionic case as shown in 



Ref.pUj. 

The main uncertainty sources for the twist-3 contribution come from the wave function 
^ p (x, kj_) and the parameter (3k- We show the contributions to the charged kaon electromag- 
netic form factor from the two groups of parameters for ^ p (x, kjj in Fig. ([6]), c.f. Eqs. fl32l 
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[331) . It is found that the uncertainty within the allowable energy region caused by these 
two groups of parameters are about 10 — 20% and 20 — 30% for the charged case and the 
neutral case respectively. Secondly, we show the uncertainty caused by the parameter (3 K in 
Fig.(J7]), where the Group 1 moments 25j are used to determine the parameters of $?(x, kj_) 
and three typical values (3k = O.SSGeV^ -1 , 0.%7GeV~ l and O^GeV -1 are adopted 4 . The 
twist-3 contribution increases with the increment of (3k, arid the uncertainty is less than 5% 
for the charged form factor, while for the neutral form factor the uncertainty changes to be 
~ 10%. 

As for the higher order corrections, we present a naive estimation of the next-to-leading or- 
der (NLO) twist- 2 contribution to the charged kaon electromagnetic form factor with the help 
of the asymptotic DA, i.e. with the renormalization scale and the factorization scale taken 
to be n 2 R = n 2 f = Q 2 , it can roughly be expressed as [39j, Q 2 F^° « (0.903GeV 2 )^ra 2 s (Q 2 ). 
Numerically the NLO correction will give about ~ 20—30% extra contribution to the charged 
kaon electromagnetic form factor. 



IV. SUMMARY AND CONCLUSION 



The kx factorization formalism provides a convenient framework and has been widely 
applied to various processes. In this paper we present a systematical study on the kaon 
electromagnetic form factors F K ± K o k°(Q 2 ) within the kx factorization formalism. In or- 
der to get a deeper understanding of the hard contributions at the energy region where 
pQCD is applicable, we have examined the transverse momentum effects, the contributions 
from the different helicity components and different twist structures of the kaon LC wave 
function. Our results show that the right power behavior of the hard contribution from 
the higher helicity components and from the higher twist structures can be obtained by 
keeping the kx dependence in the hard scattering amplitude. The full estimation of the 
power suppressed contributions to the kaon electromagnetic form factors Q 2 Fk+(Q 2 ) and 
Q 2 F K o(Q 2 ) is shown in Fig.®. The kx dependence in LC wave function affects the hard 
and soft contributions substantially and the power-suppressed terms (twist-3 and higher 
helicity components) make an important contribution below Q 2 ~ several GeV 2 although 

4 When varying af(lGeV) and a%(lGeV) within the region of [0.03,0.07] and [0.05,0.15] respectively, the 
value of Pk shall vary within the region of [0.856, O^GlGeV^ 1 . 
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FIG. 8: The combined hard contribution for the kaon electromagnetic form factors Q 2 Fx+{Q 2 ) 
(Left diagram) and Q 2 F K o(Q 2 ) (Right diagram). The solid line stands for the combined hard 
contribution, the dotted line, the dashed line and the dash-dot line are for LO twist-2, twist-3 and 
NLO twist-2 contributions respectively. 

they drop fast as Q 2 increasing. 

The parameters of the proposed model wave function can be fixed by the first two mo- 
ments of its distribution amplitude and the normalization condition. In this paper we have 
taken the first two moments af (lGeV) = 0.05 ± 0.02 and af (lGeV) = 0.10 ± 0.05. It is 
found that the uncertainty of the kaon electromagnetic form factor, which is caused by vary- 
ing values within the above range, is rather small. It is also found that the power-suppressed 
twist-3 contribution makes an important contribution at Q 2 ~ several GeV 2 and drops fast 
as Q 2 increasing. A naive estimation gives the NLO correction about ~ 20 — 30% extra 
contribution to the charged kaon form factor. 

The relativistic effect due to the Wigner rotation have also been applied to calculation 
the kaon electromagnetic form factor. Consequently there are higher-helicity (Ai + A2 = 
±1) components in the spin space wave function besides the usual-helicity (Ai + A2 = 0) 
components. It is shown that the higher helicity components have the same importance 
as that of the usual helicity components for the soft energy region, e.g. the probability 
of finding the valance states in the charged kaon includes two parts: (-Pig 1+A2_ °' ) = 0.562) 
for the usual helicity components and (Pig 1+X2 ~ ±1 * > = 0.339) for the higher helicity states 
for af{lGeV) = 0.05 and af (lGeV) = 0.115. By taking af {lGeV) = 0.05 ± 0.02 and 
af (lGeV) = 0.10 ± 0.05, we obtain the uncertainty of the probabilities P uS = 0.9011q;oio, 
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It is found that the hard-scattering amplitude for the higher-helicity components is of order 
1/Q A which is the next to leading contribution compared with the contribution coming 
from the ordinary helicity component, but it can give sizable contributions to the kaon 
electromagnetic form factor, especially for the twist-2 shown in Fig.®. 

Acknowledgements 

This work was supported in part by the Natural Science Foundation of China (NSFC) 
and by the grant from the Chinese Academy of Engineering Physics under the grant 
numbers: 2008T0401 and 2008T0402. 



[1] T. Huang and X.G. Wu, Phys.Rev. D70, 093013(2004). 

[2] T. Huang, X.G. Wu and X.H. Wu, Phys.Rev. D70, 053007(2004). 

[3] X.G. Wu and T. Huang, Int.J.Mod.Phys. A21, 901(2006). 

[4] H.M. Choi and C.R. Ji, Phys.Rev. D59, 034001. 

[5] BW. Xiao, X. Qian and B.Q. Ma, Eur.Phys.J. A15, 523(2002). 

[6] S.D. Drell and T.M. Yan, Phys.Rev. Lett. 24, 181(1970). 

[7] S.J. Brodsky, SLAC-PUB-2447, 1979; G.P Lapage and S.J. Brodsky, Phys.Lett. B87, 
359(1979). 

[8] G.P. Lepage, S.J. Brodsky, T. Huang and P.B. Mackenzie, in Particles and Fields-2, Pro- 
ceedings of the Banff Summer Institute, Banff, Alberta, 1981, edited by A.Z. Capri and A.N. 
Kamal (Plenum, New York, 1983), P83. 
[9] G.P. Lepage and S.J. Brodsky, Phys.Rev.D22, 2157(1980); ibid. 24, 1808(1981). 
[10] M. Beneke and Th. Feldmann, Nucl.Phys. B592, 3(2001); Z.T. Wei and M.Z. Yang, Nucl.Phys. 

B642, 263(2002). 
[11] W.M. Yao, et al. (Particle Data Group), J.Phys.G 33, 1(2006). 

[12] S.J. Brodsky, T. Huang and G.P. Lepage, in Particles and Fields-2, Proceedings of the Banff 
Summer Institute, Banff, Alberta, 1981, edited by A.Z. Capri and A.N. Kamal (Plenum, 
New York, 1983), P143; T. Huang, in Proceedings of XXth International Conference on High 



23 



Energy Physics, Madison, Wisconsin, 1980, edited by L.Durand and L.G. Pondrom, AIP 

Conf.Proc.No. 69(AIP, New York, 1981), plOOO. 
[13] H.N. Li and G. Sterman, Nucl.Phys. B325, 129(1992); T. Huang and Q.X. Shen, Z.Phys. 

C50, 139(1991); F.G. Cao, T. Huang and C.W. Luo, Phys.Rev. D52, 5358(1995). 
[14] T. Huang, B.Q. Ma and Q.X. Shen, Phys.Rev.D49, 1490(1994). 

[15] G. Sterman, Phys.Lett. B179, 281(1986); Nucl.Phys. B281, 310(1987); S. Catani and L. 

Trentadue, Nucl.Phys. B327, 323(1989); Nucl.Phys. B353, 183(1991). 
[16] H.N. Li and H.L. Yu, Phys.Rev.Lett. 74, 4388(1995); Phys.Lett. B353, 301(1995); Phys.Rev. 

D 53, 2480(1996). 

[17] T. Kurimoto, H.N. Li and A.I. Sanda, Phys.Rev. D65, 014007(2002); H.N. Li, Phys.Rev.D 

66, 094010(2002). 
[18] J. Botts and G. Sterman, Nucl.Phys. B225, 62(1989). 
[19] E. Wigner, Ann.Math. 40, 149(1939). 

[20] X.G. Wu, T. Huang and ZY. Fang, Eur.Phys.J. C52, 561(2007). 

[21] X.G. Wu, T. Huang and ZY. Fang, arXiv: 0712.0237, to be published in Phys.Rev. D. 
[22] H.M. Choi and C.R. Ji, Phys.Rev. D74, 093010(2006). 
[23] X.H. Guo and T. Huang, Phys.Rev. D43, 2931(1991). 

[24] C.R. Ji, P.L. Chung and S.R. Cotanch, Phys.Rev. D45, 4214(1992); H.M. Choi and C.R. Ji, 

Phys.Rev. D75, 034019(2007). 
[25] P. Ball, V.M. Braun and A. Lenz, J.High Energy Phys. 0605, 004(2006). 
[26] A. Khodjamirian, Th. Mannel and M. Melcher, Phys.Rev. D70, 094002(2004). 
[27] P. Ball and M. Boglione, Phys.Rev. D68, 094006(2003). 
[28] V.M. Braun and A. Lenz, Phys.Rev. D70, 074020(2004). 

[29] P. Ball and R. Zwicky, JHEP 0602, 034(2006); V.M. Braun and A. Lenz, Phys.Rev. D70, 
074020. 

[30] V.M. Braun et al, Phys.Rev. D74, 074501(2006). 

[31] PA. Boyle et al, Phys.Lett. B641, 67(2006); |hep-lat/0610025] 

[32] M.A. Donnellan et al, arXiv: 0710.0869. 



[33] V.M. Braun, etal, QCDSF/UKQCD collaboration, |hep-lat/0610055| Phys.Rev. D74, 
074501(2006). 

[34] Seung-il Nam and Hyun-Chul Kim, Phys.Rev. D74, 076005(2006). 

24 



[35] V.L. Chernyak and A.R.Zhitnitsky, Nucl.Phys.B201, 492 (1982); Phys.Rep.112, 173(1984); 

Nucl.Phys.B246, 52(1984). 
[36] Seung-il Nam and Hyun-Chul Kim, Phys.Rev. D74, 096007(2006). 

[37] S.R. Amendolia eta/., Phys.Lett.B178, 435(1986); E.B. Dally eta/., Phys.Rev.Lett. 45, 
232(1980). 

[38] G. Dillon and G. Morpurgo, Europhys.Lett. 54, 35(2001). 

[39] B.Melic, B. Nizic and K. Passek, Phys.Rev. D60, 074004(1999); and reference therein. 



25 



